Conductance suppression in normal metal-superconductor mesoscopic strucutres 
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o\ : 

Using a scattering matrix approach and quasiclassical Green's function technique, we calculate 
» I . the conductance of the S/N system (see Fig.l). We establish that the difference between the su- 

perconducting and normal state conductance (5G — G s — G n ) is negative for large S/N interface 
resistances (Rs/n) an d changes sign with decreasing Rs/n- The comparision of the results obtained 
with experimental data is carried out. 

Pacs numbers: 74.25.fy, 73.23.-b, 72.10.-d 72.10Bg, 73.40Gk, 74.50.+r 

I. Introduction 

Recent studies of transport properties of mesoscopic, normal metal-superconductor (N/S) structures (see reviews 
jl],^| ), have revealed a number of new physical phenomena. Examples include the subgap conductance measured of 
Vh ■ SIN (superconductor-insulator-normal metal) junctions |3[j^] , oscillations in the magneto- conductance of N/S systems 
with normal or superconducting loops |7|-|l2|,and the nonmonotonic dependence of the conductance on temperature 
and voltage Although the majority of experimental results have been successfully explained, there are some 

which remain anomalous. In particular the increase in resistance of diffusive N/S systems in a certain temperature 
range below T c [pd|-tj6| has remained unexplained for a number of years. An early theoretical prediction theit 
superconductivity induced conductance suppression is a generic feature of N/S nanostructures was followed by quan- 
titative theories of this effect in the ballistic and Anderson localisied regions H) as well as in resonant structures 
However a quantitative theory in the diffusive region has remained elusive. Several authors have suggested 
possible explanations of this puzzling phenomenon. In the simplest, |l4[ , the resistance change <5i? s = R s — R n is 
determined by a change in the interface resistance Rb which is larger in the superconducting state than in the normal 
state. Another possibility is presented in Jl^,|l^l where a two dimensional, multiprobe geometry was considered. The 
I ' current I a passes through two contacts on one side of the normal film contacting a superconductor and the voltage Vj, 
. is measured between two probes located on the opposite side of the normal film. The authors of Ref Jlj],[l8| showed 
that the quantity R a b = Vb/Ia may exhibit an increase below T cs compared with its normal state value. In this 
geometry the spatial distribution of the current is nonuniform. However in some experiments, the geometry is almost 
one dimensional with the current distribution across the width of the normal film almost uniform. Therefore this 
f~>) mechanism may not be responsible for all the experimental observations of enhanced resistance. 

In this paper we suggest an alternative mechanism which determines the change in resistance SR S (or the conduc- 
CTN . tance SG S « —SRs/R^) of the structure shown in Fig. 1. We will show that SR S may be positive if the interface 
resistance Rb is large enough compared with the resistance of the metallic film in the normal state R n . The variation 
5R S is determined by two factors: a variation of the shunting interface resistance SRb, leading to a positive change in 
resistance and a variation of the normal film resistance due to a condensate induced by the proximity effect. In view 
of these conflicting effects it is not obvious what sign SR S will adopt for any given parameters of the system. In what 
follows we use two methods to study the change in resistance, namely an analytical quasiclassical technique and a 
numerical scattering approach. The scattering approach Q complements the quasiclassical method, and enables us 
to probe areasof parameter space which lie outside the region of validity of the latter. 

II. Quasiclassical Theory 

Consider a diffusive regime where the mean free path is shorter than any other characteristic length in the system 
(except the Fermi wavelength). Such a case is realized in most experiments performed on metallic films or on doped 
semiconductors. For diffusive S/N mesoscopic structures, equations for the quasiclassical Green's functions were 
derived many years ago and are presented (in the most convenient form suitable for the present analysis ) in Larkin 
and Ovchinnikov's paper ]l9] . These equations must be supplemented by boundary conditions at the S/N interface 
derived by Zaitsev |^(| (see als o [pl|j22| ]) and have been used extensivly for the theoretical study of transport properties 
of S/N mesoscopic structures p3|-|33|]. 

In this paper we shall assume that the proximity effect is weak, i.e. the amplitude of the condensate induced in 
the normal film is small. We will show that this is true for structures where the S/N interface resistance (Rb) , in 
the normal state, exceeds the resistance of the normal film. When this condition is satisfied, the condensate functions 
obey the linearized Usadel equation and the distribution function obeys the equation (see ]2^,|30|), 



1 



L 2 d x ((l - m)d x f) = fg b G b d(x e (S/N)), (1) 

where the function $(x) is equal to 1 in the S/N region and zero otherwise, m = \Tr(F R - F A f and F R ^ is the 
retarded (advanced) Green's function, g b = pL 2 / R b nd = g b \{LjL\) , g b \ is the ratio of the normal film resistance to 
the S/N resistance, R b o is the S/N imterface resistance per unit area in the normal state, p and d are the specific 
resistivity and the thickness of the normal film. The function G b {x) determines the local normalized conductance of 
the S/N interface in the superconducting state, 

G b (x) = v s v n + ±Tr(F R + F A )(F S R + F A ), (2) 

where the density of states (DOS) in the superconductor v s — Re(e + iT)/^/((e + iT) 2 — A 2 ) and v n is the DOS 
in the normal film (for simplicity we assume v n = 1 ,i.e. the S and N metals are regarded as identical apart from 
the critical temperature, we also assume that T cn = 0).The first term in Eq. (^) describes the contribution of the 
quasiparticle current to the conductance (if Y = it differs from zero only at energies |e| > A). The second term 
is due to Andreev reflection and discribes a conversion of the low-energy quasiparticle current into the condensate 
current (if T = 0, the current is not zero for |e| < A). The condensate functions F^ A ^ in the superconductor are 
assumed undisturbed by the proximity effect (this is true provided that A » 75 = K 2 D , where K 2 = (Rbada)^ 1 ), 
and they are equal to, 

F s R M=iT y F s R W, (3) 

where F^ A ^ — A/^/((e + iT) 2 — A 2 ), T is the damping rate in the excitation spectrum of the superconductor. 
Assuming that the right-hand side of Eq. ([!]) is a small perturbation we easily find a solution for, 

t _ / {JA X + Jo dxi(m + mb)]) , < x < L\ , . 

J I (J(x-L 1 )+/(L 1 )), L l<x <L W 

where J\ and J are the energy dependent intergration constants. The current I through the system is expressed in 
terms of J ( see @|§), 

I = (ad/2e) [ deJ{e). (5) 



The function mb in Eq.(Q) is given by, 

nib = (gb/L 2 ) dx 1 x 1 G b (x 1 ). (6) 
Jo 

In the reservoirs the distribution function has the equilibrium form, 

/ (e, L)=F V = [tanh (e + eV) - tanh (e - eV) 0] /2, (7) 

where = (2Tkg)^ 1 . By matching the functions and their derivatives at x = L\ and using Eq. ([?]), we find for the 
"partial current" J(e), 

J(e) = (F v /L)[l- < m > +h(m bl - < m b >i)], (8) 

The angle brackets mean a spatial averaging over the regions (0,L) and (0, Li);mb\ = mb(Li)0i = L\/L. With 
the aid of Eqs. (|) and (|), we find the normalized difference between the differential conductances of the system in 
the superconducting state and the normal state, 

SS = Gs ~ ° n = - f de0F' v [< m > +h(< m b - m n h > x -(m M - m? x ))]/[l + g b l\l% (9) 

G n Jo 

Here G s<n = (dl /dV) Si . n is the differential conductance below and above T cs ; is the function m b (x) in the normal 
state: m£ = g b x 2 /2L 2 , F^ = 8F v /d(eV0). 

Let us discuss the physical meaning of the different terms in Eq. (|^) . The first term gives a positive contribution 
to SS (< m > is negative). This term arises from the renormalization of the normal film conductance caused by the 
induced condensate. This has been calculated in several papers |30|-|33|], where it was established that this term has 
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a nonmonotonic voltage and temperature dependence decreasing to zero at [eV, T] = and [eV, T] 3> £l (cl = D/L 2 
is the Thouless energy). The second term in Eq. (^|) determines the change in conductance due to different values 
of the S/N interface resistance in the normal and superconducting states. The contribution of this term to SS is 
negative because the S/N interface resistance in the superconducting state is larger than the normal state interface 
resistance (as long as the barrier transparency is not too high). Let us estimate the magnitudes of these terms. If 
g b is small, the second term in Eq. (|^) (the subgap conductance) is small compared with the first term. Therefore 
at low temperatures the contribution caused by the second term in Eq. (^) is related to a change in the DOS of the 
superconductor. This yields, 

SS dos » -g b ll/3. (10) 

As we shall see, the amplitude of the condensate functions F R ( A ' induced in the normal film by the proximity 
effect is of the order g b l\ , i.e. of the order of the ratio of the normal film and the S/N interface resistances. The 
characteristic energy of decay of F R ^(e) is the Thouless energy e^. Thus the contribution to SS from the proximity 
effect (the first term in Eq. (]|)) is (if e L < T < A and h < 1), 

SS pr « gill (11) 

Comparing Eqs. ( |Io|) and (|ll|), we see that the sign of SS changes from negative to positive as g b increases and 
5S m in « —if is reached when g b ra l\. 

In order to find SS, we need to calculate F R ( A \ As noted above, in the limit of a weak proximity effect the function 
pR( A ) obeys the linearized Usadel equation which may be presented in the form (see for example j2j],|3(J ) , 

d xx F R ^ - (k R ^) 2 F R ^ = -{g b /L 2 )F R ^d(x g (S/N)), (12) 

where (k^'^) 2 = (t2«c + j)/D, 7 is the depairing rate in the normal film. The solution to Eq. ( p"2| ) satisfying the 
boundary conditions F R ^ (±L) = is the function, 

F = (g b /9 2 )F s ( S 1 7 ZTu ikX \\{ ° T <X<L l (13) 

Here 9 = kL, C2 = cosh^)/ cosh(6>),si = sinh(6*i)/ cosh(0), #1.2 = kL\ i, L2 = L — L%. For convenience we have 
dropped the indices R(A). One can see from Eq. ( [l3| ) that at characteristic energies e k t j, « A the amplitude of 
pR( A ) i s f the order g b l\. It is worth noting that if l\ <C 1 (this condition is satisfied in most experiments), then the 
magnitude g b l\ corresponding to the actual g b is of the order Z^,i.e. as l\ is small, the proximity effect is small. 

From Eq.©, we can find the quantites < m >,m b \ and < m b \ >. Substituting them into Eq. (|9|), we obtain the 
variation of the normalized conductance 5S as a function of tcmperature,f/i, l\ etc (see Appendix). In the general 
case the expression for SS has a rather complicated form, but may be simplified drastically by taking the zero-bias, 
zero-temperature limit, in which case 5S becomes, 

SS =g b l 3 1 [g b h(l-h^-)-\}. (14) 
15 3 

In this case the contribution due to the variation of the normal region conductance goes to zero (the term < m > 
in Eq.(^)). The variation SS is caused by a change in the S/N interface conductance G b ; the first term in Eq.(||) 
gives a negative contribution, and the second term (i.e. the subgap conductance) gives a positive contribution. This 
expression also changes sign at g b l\ w 1/3. However in this case the condensate amplitude is not small, and strictly 
speaking, Eq. ( [Ti]) is not valid when g b l\ is of order 1. Nevertheless, we show later using numerical calculations that 
this conclusion regarding the change in sign of SS remains valid in the zero-bias, zero temperature limit. Fig. 2 shows 
the dependence 5S(g b ) for l\ = 0.2 and l\ — 0.4. We see that in accordance with qualitative speculations given above, 
SS is negative at small g b , reaching a minimum with increasing g b and then changing sign. The magnitude of SS m i n 
decreases with decreasing l\ and depends on temperature in a complicated nonmonotonic way (see Fig. 3). 

In Fig. 4 we show the dependence of SS on 7 (7 may increase by applying an external magnetic field, 7 ~ 
H 2 ), the effect of a negative SS becomes more pronounced. This behaviour is quite clear from a physical point of 
view. The applied magnetic field suppresses the proximity effect, but affects the DOS only weakly. Therefore the 
relative contribution SSdos increases with increasing magnetic field. A similar situation takes place in fluctuation 
paraconductivity in layered superconductors where an increase in the resistance due to superconducting fluctuations 
is enchanced by the magnetic field Q . 
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In Fig. 5 we plot the voltage dependence of SS. One can see that two maxima exist in this dependence; one is close 
to zero bias and another one located at eV w A. A similar voltage dependence of the phase coherent conductance 
has been observed in the recent work of |3f|. We note that although the effect of the conductance decrease is small 
(< 5%), in dimensional units the conductance decrease SG may be much larger than the quantum conductance 2e 2 /h 
( 8G 3> 2e 2 /h). Only in this limit can we use quasiclassical theory. 

III. Numerical simulations 

In this section we use the scattering approach reviewed in Q to determine SS ,for a tight binding lattice with the 
geometry of Fig. 1. In the linear-response limit, at zero temperature, the conductance of a phase-coherent structure 



may be calculated from the fundamental current voltage relationship 37 



h \ _ ( an a u \ ( V\ - v 

h) I fl 21 022/ \V 2 -V 



(15) 



where, at finite temperatures, 



f an 




2e 2 r 


\CL2l 


a 2 2 j 





f]F( _dm ) fN+(E)-Ro(E)+Ra(E) T' a {E)-T'{E) \ 

[ dE ] \ T a (E)-T (E) N+(E)-R> (E) + R> a (E))> 

Eq. ( |l5| ) relates the current Ii from a normal reservoir i to the voltage differences (vj — v), where v = fi/e (fi is the 
chemical potential of the superconductor). The ay's are linear combinations of normal (To, Rq) and Andreev (T a , R a ) 
scattering coefficients. The primes on the coefficients refer to quasiparticles originating from the right-hand reservoir, 
whilst the coefficients without primes refer to particles from the left reservoir. Setting I\ = I = —1% and solving Eq. 
( |l5| ) the two probe conductance is (see @), 

q _ aii<Z22 — ai2<221 

an + a-22 + ai2 + fl2l 

As noted in |59| in the presence of disorder, the various transmission and reflection coeffcients can be computed 
by solving the Bogoliubov - de Gennes equation on a tight-binding lattice of sites, each labelled by an index i and 
possessing a particle (hole) degree of freedom ip(i) (<p(i)) (ip(ip) is the particle (hole) wavefunction) . In the presence 
of local s-wave pairing described by a superconducting order parameter Aj, this takes the form, 

Elpi = Ci^i -Y,6 T (^i+S + i>i-6) + AiVi (18) 
Eipi= -em +J2s T (<Pi+6 + Vi-s) + ^i^i- 

In what follows, in the normal diffusive region, the on-site energy is chosen to be a random number, uniformly 
distributed over the interval eo — 1 to eo + 1 , whereas in the clean N-regions ej = eo. In the S-region, the order 
parameter is set to a constant, A,; = Ao, while in all other regions, Aj = 0. The nearest neighbour hopping element r 
merely fixes the energy scale (i.e. the band- width), whereas eo determines the band- filling. In what follows we choose 
t = 1. By numerically solving for the scattering matrix of Eq. (fl8| ), exact results for the dc conductance can be 
obtained p8| , p7y3q |. In the zero bias, zero temperature limit, Eq. ( |l7|) is greatly simplified and reduces to, 

G = T +T a + 2{RaR a- T - T a) _ (1Q) 
Ra + R a + T a + T a 

For the structure shown in figure 1 , with a superconductor of length 2Li , and a barrier resistance R, evaluation 
of this expression yields results for < G n >,< G s > and < SG > shown in table 1. In each case, the normal 
diffusive region is 40 sites wide and 64 sites long. The superconductor is of width 20 sites with Ao = 0.1 (Ao = 0) 
in superconducting (normal) state. Results are obtained by averaging over 100 disorder realizations, yielding an 
estimated error in the mean values of approximately 0.04. the first row of the table shows results for Li = 30, R = 2 
and demonstrates that a negative SG can indeed occur, and comparison with the R = 0.5, shows that lowering the 
interface resistance causes SG to change sign. As discussed previously, this result was expected although could not be 
proved using quasiclassical theory. Also, as discussed previously, when Li is decreased (e.g. to Li = 20 with R = 2) 
the table shows that SG changes sign and becomes positive. Finally, to examine the effect of a magnetic field, the 
forth row of the table shows results with a magnetic field applied to the normal region (corresponding to 0.8 flux 
quanta through the whole structure). This demonstrates that the introduction of a magnetic field causes a negative 
SG to become more negative in agreement with the quasiclassical approach but in conflict with experimental evidence. 
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3.79 


3.47 


-0.33 



Finally we note that at a finite temperature (ksT = cl) , where the full integral of Eq.(|lq) needs to be evaluated, 
we find for the structure of row 1 in the table, < G n >= 3.68, < G s >— 3.56,< SG >= —0.12, which confirms 
the prediction made using the quasiclassical approach, that the onset of superconductivity causes a drop in the 
conductance of the structure, even at finite tempretures. 

IV. Discussion 

We have demonstrated that superconductivity-induced conductance suppression is an inherent property of the 
structure of figure 1. The suppression of the conductance at temperatures below T cs , is < 5% and it is enhanced by 
the application of a magnetic field. In the experiment of jl6| a stronger effect (10 — 20%) is observed, which decreases 
when a rather weak magnetic field is applied. This suggests that the magnitude of the effect and the field dependence 
are geometry-dependent. For the future, it would be of interest to confirm this experimentally by measuring the 
conductance of S/N structures of the type shown in Fig. 1 with different ratios of the normal channel and S/N 
interface resistances. 
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VI. Appendix 

Using Eq. ( |l2| ) we easily find the expression for the quantites in Eq. (^) . We have, 

< m >= ~\l\F s \ 2 < \b\ 2 > +Re(< b 2 > F 2 )}, (20) 



where F s is defined in Eq.[5[ 

g» (I h i M 2 , sinh(2^) Bjgh(2jig0 OD , sinhgig) H 2 , sinh(2/ 2 fl') gmg^Q 
< \b\ >= W (h[l + —(^qF- + 2h9 „ ) - 2i?e(c 2 — + ; 2 — (^-^ ZM*-^ 



(21) 



9b n n | 4 smh(2h9) sinh(lifl) sf sinh(2l 2 fl) 



,2r ^ ~ 1 , ^z 1 , sinh(0Zi) cosh(^i) - 1 
TOfci - m bl = g b li[—j— + —Re-^(- - c 2 ( — — )], (23) 



<m b>1 -< m? > 1= + - ^ (cosh^M + 1 - (24) 

^ = ^gjU = ^S] 1 .« - *M = \F S \ 2 - Fl 9' = Re{9},9" = Im[9}. 
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FIG. 1. The structure considered. 
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FIG. 3. Results showing dependence of 5S on temperature for gb = 1.5, 0.7, T — 0.1, h — 0.4, A = 10, 7 = and V = 0. 
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FIG. 4. Results showing dependence of SS on g b ,for T = 0.1, h = 0.4, A = 10, 7 = 0, 0.4, 0.8, and V = 0. 




FIG. 5. Results showing dependence of SS on V,for T = 0.1, h = 0.4, A = 10, 7 = 0, gb = 1.5, a = 1.0 . 
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